Analyzing temporal and spatial forest carbon storage using Google Plus Code: a case study of Zijin Mountain National Forest Park, China

Background It is always a challenging job to compare forest resources as there is not a standardized spatial unit with location information. Google Plus Code, the newest alphanumeric geocoding system, uses 20 specifically selected letters and numbers to assign a unique global ID to every cell at different levels of a hierarchical grid system which is established based on latitude and longitude. It can be used as a standardized, unique global geospatial unit to segment, locate, quantitate, evaluate, and compare natural resources with area, boundary, and location information embedded. Results For this proof-of-concept case study, forest inventory data from 1987, 2002, and 2019 for the Zijin Mountain National Forest Park in Jiangsu Province, China was analyzed based on Google Plus Code grid/cell. This enabled the quantification of carbon storage at each cell allowing for the comparison of estimated carbon storage at same or different locations over time. Conclusions This methodology is used to quantify the impacts of changing forest conditions and forest management activities on carbon storage with high spatial accuracy through the 32-year study period. Furthermore, this technique could be used for providing technical support and validation of carbon credit quantification and management. Supplementary Information The online version contains supplementary material available at 10.1186/s13021-024-00258-0.


Background
Sustainable forest management at the global, national, and regional scales requires timely and accurate forest inventory data that are georeferenced using a standardized spatial unit with the appropriate location information embedded.This need is particularly acute when it comes to measuring and estimating forest carbon storage for climate change mitigation strategies [1][2][3][4].Carbon storage is commonly estimated based on forest biomass volumes or growth rates being multiplied by the appropriate carbon index [5][6][7].Beyond total carbon storage estimates, the spatial distribution of the carbon stored within a forest is necessary to manage forests for Carbon Balance and Management optimum carbon storage [1,6,8].From these spatial distribution patterns, a cartographical presentation of forest carbon storage can be useful to provide stakeholders and decision makers better tools for making the most appropriate policy decisions or economic investments.However, these cartographic analyses require widely available and accurate georeferenced spatial units.
Commonly used spatial units like street addresses are limited for defining stand conditions in forestry applications and latitude/longitude are applicable only to specific spatial points [9][10][11].Techniques like the Kriging geostatistical method have been successfully used to develop forest carbon storage maps [8].In addition, the spatial error model (SEM) has been employed to relate the distribution of forest carbon storage to both forest stand parameters and topographic characteristics [12].Geospatial regression has also been employed as a tool for relating the spatial distribution of biomass estimates and carbon storage in forests [6].While Kriging, SEM, and geospatial regression have been applied for cartographic representation of forest carbon storage, these previous carbon storage studies usually analyze a large area of forest all together without breaking the study area into comparable and geolocatable spatial units, thus they are not helpful for on the ground carbon storage monitoring and management.
A novel technique has recently become available for a wide variety of geospatial applications, the use of Google Plus Code as an alphanumeric geocoding hierarchical grid system.This system provides exact spatial locations and unique global identification over wide ranging areas.When comparing with traditional unregular plot or subcompartment boundaries, Plus Code cells provides foresters a more convenient and accurate geolocating way for on the ground monitoring and management.However, the applicability of this grid system for forest carbon storage has yet to be developed and demonstrated for temporal and spatial forest dynamics.Therefore, this research project was initiated utilizing forest inventory data from 1987 to 2019 in one of China's most iconic scenic natural areas, Zijin Mountain National Forest Park.These data were employed to determine georeferenced carbon storage utilizing Google Plus Code, with the purpose of providing researchers and decision makers with an effective method for long-term forest carbon stock management.

Study area
Zijin Mountain National Forest Park is located in Nanjing, Jiangsu Province, China.The area's latitude and longitude ranges from 32°01'57"N to 32°18'15"N and from 118°48'24"E to 118°53'04"E, respectively.The highest elevation in the park is 448.9 m and it has an area of

Forest inventory field data collection
The forest in this research is where there is a canopy within the park and forest management inventory were conducted in 1987 and 2002.These historic data were obtained from the Zijin Mountain National Forest Park administration database.From May to December in 2019, a third forest inventory was conducted.In total, there were 777 sub-compartments or stands delineated within the study area (Figs. 1 and 2).Sampling or measurement plots were 666.67 m 2 (25.82 m x 25.82 m) with approximately one plot per three ha, based on compartment size.At each measurement plot, aspect, slope, soil type, soil depth, species, age, age class, diameter at breast height (DBH), height, and canopy closure were recorded.Forest condition and tree volume were then estimated for each sub-compartment based on these measurements.

Estimation of forest biomass and forest carbon storage
For estimation of forest biomass, the method described by Fang et al. [4] was used with the biomass conversion continuous function with coefficients being described by: Where B is the forest biomass (t, metric ton), V is the forest volume (m 3 ), a and b are conversion coefficients which change depending on forest type or tree species.Measurements of litter, understory, and grass biomass were not available so only above-ground biomass was estimated in this study.The a and b values are adopted from researches done by Zhang et al. [13], Zhao et al. [14] and Fang et al. [5] and are summarized in Table 1.Forest carbon content (q i , Table 1) was obtained from the Guide for Carbon Accounting and Monitoring in Afforestation Projects which utilized the Intergovernmental Panel on Climate Change (IPCC) reference values and bamboo forest carbon content of 0.50 [15,16].The formula used to calculate forest carbon storage (C) is: Forest carbon density is the forest carbon storage per unit area and was calculated as: ρ is forest carbon density (tC/ha), C is carbon storage (tC) and S is the corresponding spatial area (ha).

Geospatial distribution of forest carbon storage using Google Plus Code
There is currently not a standard geospatial unit with embedded location information that is widely utilized for forest management.In recent years, alphanumeric geocoding systems, such as What3Words, Mapcode, and Placekey, were developed to locate any point on Earth.The newly available Google Plus Code has a similar Although none of these geocoding methods were designed for forestry applications directly, the hierarchical grid system of Plus Code enables foresters to have comparable cells for the entire world.Google Plus Code uses latitude and longitude to start the level zero grid (20° by 20°) and each cell has a two-digit (letter and number) ID.Only 20 specifically selected letters and numbers are used to prevent from confusion that can occur between certain alphanumeric pairings, such as lower-case L and the number one.Each 20° by 20° macro-cell is then subdivided into level one sub-grids (1° by 1°).At this level, each cell will have a unique four-digit ID (two digits from level zero and two digits for level one itself ).This subdivision is further sub-divided until a level five (12 digits ID) is determined.A different structure is used for level six + divisions (Fig. 1).Thus, regardless of the level, every hierarchical cell on Earth has a unique global ID assigned.
In other words, Google Plus Code established a reference system with unique ID for every cell in a hierarchical grid system.Google (google.com/maps or plus.codes/map) can provide location and direction to any specific cell directly with its ID.For example, when searching "8P4W3VQC+", Google Map will show the central point of the cell and plus.codes will show the cell boundary on the map.Along a similar latitude, the cells of this welldefined grid system have same area The forest carbon storage CV index (Coefficient of Variation) can be used to represent the relative stability of a cell's forest carbon storage.The larger the CV index, the less stable the cell's forest carbon storage is.On the contrary, the smaller CV index means a more stable forest carbon storage.CV index can be calculated with the equation below: C i is the forest carbon storage at year i and C is the aver- age forest carbon storage for each cell from 1987 to 2019.Level 2, 3, and 4 grid cells have an area of 2,630, 6.6 and 0.016 ha, respectively.For carbon storage and forest management, level 3 cells of Plus Code were used as the location and unit boundary for this study.There are 588 cells and each cell has an area of 66,123 m 2 (237 m × 279 m).At this level, the ID for each cell has 8 digits, such as 8P4W3VQC (Fig. 2).For every forest inventory survey year, the carbon storage of each cell was estimated with ArcGIS 10.5® based on the methodology described in this section.Proportional method was applied for cells with multiple sub-compartments.

Forest carbon storage change
Forest area in the park (including all tree species and bamboo forest) increased from 1791.5 ha in 1987 to 2,311.1 ha in 2019 (Fig. 3).This increase in forest area is noteworthy particularly given the increasing number of visitors per year to the park and the rapid urbanization and population growth in the Nanjing metropolitan area over that time period.The rate of increase in forest area from 2002 to 2019 (207 ha) is lower than that from 1987 to 2002 (312 ha) due to the higher initial reforestation efforts that occurred in the park in the late 20th century.Mixed conifer and broadleaf forest area increases accounted for most of this change.Forest carbon storage also increased for the study area (Fig. 2).However, the rate of forest carbon storage from 2002 to 2019 is higher than that from 1987 to 2002 due to the forest type and age structure.For the 32-year study period, the net increase in forest carbon storage was 43,740.7 tC.The average annual increase in forest carbon storage was 1,197.6 tC/  2).For example, broadleaf forests accounted for just over half (51-54%) of the carbon stored per ha compared to conifer which accounted for approximately a quarter (24-28%).Bamboo remained relatively constant throughout the study period both in terms of area (2%) and the percentage of the carbon stored per ha (11-14%).The mixed conifer/broadleaf forest type increased the most in terms of forest area, from 8% (126 ha) of the park in 1987 to 31% (443 ha) by 2019.This mostly occurred as stands of pure Masson pine forests were converted into mixed forest types, with Masson pine area decreasing from 444 ha in 1987 to 32 ha in 2019 due to insect mortality.However, the mixed conifer/broadleaf forest type accounted for the lowest proportion of carbon storage, ranging from only 8-9% of the tons of carbon per ha from 1987 to 2019.This lower percentage of the per ha carbon storage can in part be explained by the younger age class of this forest cohort, as noted below.
Across all forest types, the mixed conifer class had the greatest increase in carbon storage per ha, increasing by 17.1 tC/ha from 37.4 tC/ha in 1987 to 54.5 tC/ha in 2019.For the broadleaves forest, oak had the greatest increase in rate of carbon storage, steadily increasing by 15.7 tC/ ha from 44.2 tC/ha in 1987 to 59.9 tC/ha 2019.The rate of carbon storage for Masson pine increased by 13.7 tC/ha and 8.6 tC/ha locust from 1987 to 2019.Formosan gum remained fairly consistent in its rate of carbon storage, decreasing slightly by 1.6 tC/ha.Carbon storage rates are related to stand age structure in addition to species, as noted below.Bamboo was the most consistent cover type in the park in terms of carbon storage per ha in part because of its more rapid maturation rate when compared with tree species.

Forest carbon storage by age cohort
As noted above, forest carbon storage is dependent on forest age structure.Age classification is based on tree species and cover type which were derived from province level forest management practices (Table 3).
Forests in this park are primarily managed for their scenic value and recreation opportunities for visitors and there was no forest harvesting in the study area during this 32-year study period.Overall, stand cohorts aged by 32 years and there was an estimated net increase in total accumulated carbon in the park with this increase in biomass (Table 4; Fig. 4).For example, young forest carbon density increased from 17.7 tC/ha to 22.5 tC/ha and 31.0 tC/ha for 1987, 2002, and 2019, respectively.This can be attributed to hardwood reforestation in the study area over the years.Similarly, in 1987, middleaged and mature forests covered 76.7% of the area and correspondingly, their carbon storage accounted for about the same proportion: 76.4% of the total storage.In 1987, over-mature forest occupied only a small area with little carbon storage (0.2%).By 2002 and 2019, the over-mature forest area increased as did the total carbon stored in these forests, accounting for 16.5% (2002) and 30.0%(2019) of total storage.These results are consistent with those reported by Chen et al. [18] for southeastern China.For the current study, the rate of accumulation by forest type for over-mature forests was not evaluated due to data limitations.That said, consistent with the findings of Chen et al. [18], as over-mature forests begin to senesce in growth rate, the rate of carbon accumulation per ha will likely continue to decrease.This may necessitate some forest harvesting and management to prevent these over-mature forests from releasing this stored carbon as they decay in the coming decades [19].

Analysis of forest carbon storage distribution based on Plus Code cells
As mentioned before, 588 Plus Code level 3 grid cells (237 m × 279 m) were used as the location and unit boundary for the study area.For every forest inventory survey year, the carbon storage of each cell was calculated and classified into 8 categories (Fig. 5).The entire area was then grouped into low storage area (less than 200 tC/cell), medium storage area (200-400 tC/cell) and high storage area (400-800 tC/cell) (Table 5).Forest carbon storage is greater for the northern cells which is more heavily forested since there is a greater proportion of scenic sites and tourist areas in the southern part.From a temporal aspect, carbon storage continued increasing, especially for the southern cells during the past 32 years.The coverage of medium storage area increased from 44.8% (1987) to 49.4% (2002) and 52.0%(2019).The coverage of high storage area increased from 0.3% (1987) to 4.5% (2002) and 8.3% (2019).In summary, the total forest carbon storage increased from 1987 to 2019, as would be expected given the forest restoration efforts and improved management practices in the park.
Plus Code cells were classified into 5 categories based on carbon storage CV index, with low stability (13.3%), low-medium stability (12.8%), medium stability (12.9%), medium-high stability (33.1%), and high stability (25.1%) (Table 6; Fig. 6).In terms of stability, 71.1% of the forest is   Forest conservation and management management measures should be used to maintain and enhance carbon storage along with the other ecosystem services that the forests provide in this park.In addition, from the geospatial aspect, Google Plus Code was found to be an effective grid system for analyzing temporal and spatial trends in carbon storage.Using this system, it was found that the southern portion of the study area had less carbon storage than the northern portion in 1987 and 2002.However, carbon storage increased by 2019.Medium and high carbon storage area coverage increased from 45.1% (1987) to 53.9% (2002) and 60.3% in 2019.The Google Plus Code analysis indicated that perhaps future restoration and forest protection efforts should be focused in this area.Pairing Google Plus Code with the CV index showed that the carbon storage stability is relatively high for entire study area.
One last thing worth mentioning is that Google Plus embedded with a unique global ID to which forest conditions can be locked to a particular spatial unit in time.This enables more efficient quantification of changes in land use and forest conditions.This case study was focused an important concern for China, carbon storage and sequestration through forest restoration in one of its iconic natural areas, Zijin Mountain National Forest Park.Google Plus Code can be a valuable tool to help foresters to understand and monitor the impact of forest management activities on carbon storage with high-resolution temporal and spatial accuracy.This can also provide technical support for carbon credit management.Additional studies are needed to assess the potential of Google Plus Code for the other aspects of ecosystem management, including but not limited to the biomass of litter, understory, grass, below-ground biomass, and water resources to provide a more comprehensive picture of the eco-function that forests provide.Future research can include but not limited to uncertainty analysis, applying this method to compare carbon storage at different locations at different scale, and studying the impact of forest management practices on carbon storage.

Fig. 2 Fig. 1
Fig. 2 Google Plus Code level 3 grid (8-digit) and the forest sub-compartments of the Zijin Mountain National Forest Park in Jiangsu Province, China

Fig. 3
Fig. 3 Change of forest area and forest carbon storage in the Zijin Mountain National Forest Park from 1987 to 2019

*
cd stands for carbon density

Fig. 6
Fig. 6 CV Index distribution of forest carbon storage in the Zijin Mountain National Forest Park from 1987 to 2019 3,008.8 ha.Approximately 76.8% (2,311.1 ha) is forested.Due to the mid-latitude location, the region has a typical northern subtropical climate with about 2,213 h of sunshine each year and an annual average temperature of 15.4℃.The rainy season is typically in June and July and average annual precipitation is 1,090.4mm, with 70% occurring during spring and summer and 30% during fall and winter.Soils are mostly acid or weak acid yellow brown Dystric Cambisols and Eutric Planisols.Evergreen broadleaf species are the dominant canopy trees in the forest with some deciduous broadleaf species and bamboo.The predominant species include Masson pine (Pinus massoniana Lamb), Formosan gum (Liquidambar formosana Hance) and various oak species (Quercus spp.).
It is worth mentioning that the open source nature of Plus Code enables it to be directly used within any Geographic Information Systems (GIS), such as ArcGIS® and QGIS.Level 3 Plus Code grid (8 digits ID) is used for the study area in this research in order to compare the spatial change of forest carbon storage.
and provide foresters a useful spatial unit to locate, analyze and compare forest resources at the global level.Cell size does vary at different latitudes due to different ground distances corresponding to same unit of longitude.For example, the latitude ranges from 32.041° to 32.091° for our study area, cells area change from 65,455 m 2 on the south side to 65,430 m 2 on the north side.The minor difference, 25 m 2 , is ignorable.A coefficient can be used if cells at significant different latitudes need to be compared.

Table 1
Biomass conversion equations and carbon content rates for forest cover types in the Zijin mountain national forest park in Jiangsu Province, China [17]s the forest biomass (t), V is the forest volume (m 3 ), B b is the bamboo biomass (t), S b is the total area of bamboo forest (hm 2 )[17]year and 1,365.5 tC/year for the period of 1987 to 2002 and 2002 to 2019, respectively.Changes in forest carbon corresponded to different forest types and varied among measurement periods (Table

Table 2
Carbon storage and carbon density of different forest types in the Zijin Mountain National Forest Park from 1987 to 2019

Table 3
Forest age classification (years) for forest cover types in the Zijin Mountain National Forest Park in Jiangsu Province, China in the category of medium stability or greater, indicating that the forest carbon storage of the entire park is relatively stable.Most of the cells with low stability cells are on the south facing slopes, where forest carbon storage increased most from 1987 to 2019.Cells with relatively small CV Index tend to be located on the north facing slope where rates of change have been lower.Thus, forest carbon storage changed less and is relatively stable.

Table 4
Carbon storage and carbon density of different age groups in the Zijin Mountain National Forest Park from 1987 to 2019